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ABSTRACT 

The modal element method is used to study acoustic scattering from ellipses, which may be 
acoustically soft (absorbing) or hard (reflecting). Because exact solutions are available, the results provide a 
benchmark for algorithm performance for scattering from airfoils and similar shapes. Numerical results for 
scattering from rigid ellipses are presented for a wide variety of eccentricities at moderate frequencies. These 
results indicate that the method is practical. 


INTRODUCTION 

Acoustic and electromagnetic scattering from airfoils is of interest to the aerospace industry in, for 
example, studies of noise abatement and identification by radar signature. The study of scattering from 
cylinders and bodies embedded in circular cylinders using the modal element method was recently presented 
by Baumeister and Kreider (1993a). The modal element method is a computationally efficient formulation in 
which both analytical and finite element solutions are combined. This study is concerned with scattering from 
ellipses using the modal element method. The elliptic geometry can be computationally efficient in many 
aeronautic applications. For example, thin airfoil shapes can be more naturally embedded in an ellipse than in 
a circle. 

The modal element method has been employed in electromagnetic and acoustic scattering and duct 
transmission problems (Baumeister and Baumeister (1994), Baumeister and Kreider (1993a), Baumeister and 
Kreider (1993b), Lee & Cendes (1987), Astley and Eversman (1981), and Chang & Mei (1976), for example). 
In this method, the acoustic pressure or electromagnetic field in and on the scattering body is represented by 
finite elements (the finite element solution), while the field is represented by an eigenfunction expansion 
away from the body (the modal solution). The two solution forms are coupled through an interface condition. 

The purpose of this paper is to present the results of the study of acoustic scattering from ellipses 
using an elliptic interface rather than a circular one. The modal element method using an elliptic interface 
can be used effectively to study scattering from airfoil shapes for two reasons. First, if the far field results 
match the exact solution for hard ellipses of high eccentricity, it can be inferred that the method can be 
applied to scattering from a similar shape, such as an airfoil. Second, airfoils can be more efficiently 
embedded in an ellipse than in a circle, which decreases the size of the finite element domain, saving 
computation time and lessening computer memory requirements. For example, a reduction of grid size by 10 
yields a savings in the global matrix size of 100. 

There have been many recent studies of scattering from ellipses and more complicated shapes using 
other analytical methods; among them are Kim (1991), Jin and Liepa (1988), and Rojas (1988). The present 
work adds to these studies by allowing more flexibility in setting absorbing material around the surface of the 
ellipse. 
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The numerical results presented here are for scattering from hard ellipses. The modal element method 
could also be used tor scattering from soft ellipses, with the modifications mentioned under Results. 


SOLUTION METHOD 

The goal here is to compute numerically the acoustic scattering of a plane wave, travelling in the +x 
direction, that strikes an ellipse. The spatial domain is divided into two subdomains, the finite element 
domain, which contains the body, and the homogeneous domain, which surrounds the body and extends to 
infinity (fig. 1). Linear triangular elements are used in the finite element domain to calculate the pressure at 
the nodes (the finite element solution). In the homogeneous domain, an eigenfunction expansion with 
unknown modal coefficients represents the acoustic pressure (the modal solution). The two solution forms are 
coupled by imposing continuity on the pressure and velocity at the elliptical interface between the two 
subdomains. This coupling results in a single matrix equation in which the eigenfunction coefficients and the 
pressures at the finite element nodes are calculated simultaneously, yielding a global representation of the 
acoustic field. 


GEOMETRICAL MODEL 

For the numerical studies done here, the scattering ellipse is taken to be rigid; in this case, it can be 
represented by a thin ring of finite elements as shown in figure 2. The elliptic co-ordinate system used is shown 
in figure 3. The variables are standard: 


jc = — coshw cosv 
2 

y = — coshu sinv (1) 


z-z 


E, = coshw, T| = cosv 


In (u,v) co-ordinates, u is the radial and v is the angular component; in (£/n) co-ordinates, £, is the radial and 
r) is the angular component. The semi-major and semi-minor axes are denoted a and b, respectively, and 
d = 2 (a 2 -b 2 )l , 2 is the interfocal distance. The major axis is oriented along the x axis. The eccentricity of the 
ellipse is defined as 


d 1 

2a~ 


( 2 ) 


GOVERNING EQUATION 

Acoustic propagation in two dimensional space can be modeled by the continuity, momentum, and 
state linearized gas dynamic equations in the absence of flow. For harmonic pressure propagation in an 
inhomogeneous material, the following dimensionless wave equation applies: 
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The harmonic time dependence e- ia)t has been factored out. In analogy with electromagnetic scattering, £ is 
the acoustic "permittivity ", p is the acoustic "permeability" and CO is the dimensionless frequency. Details of 
these acoustic properties can be found in Baumeister and Dahl (1987). The wave number is 


k = ^ 


co 2 |i£. 


( 4 ) 


At the interface between the finite element region and the analytic region, continuity is imposed on 
the pressure and velocity. The radiation boundary condition at infinity is automatically satisfied by the 
eigenfunction expansion introduced in the next section. 


ANALYTIC SOLUTION 

In the homogeneous domain, an exact eigenfunction expansion can be derived from equation (3) by 
separation of variables. In elliptic co-ordinates, the eigenfunctions are the Mathieu functions. A plane wave 
with incident angle <t>o that strikes the rigid ellipse generates a pressure field written as (Bowman, et al., 

(1969), pg. 146, equation (3.76)) 
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p l is the incident plane wave, p s is the scattered field, c = kd/2 and 0 O is measured counterclockwise from the 
+x axis, represents the surface of the ellipse, while the observation point (£,T|) is exterior to the ellipse. The 


even and odd radial Mathieu functions of the first 
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respectively, to the Bessel functions J m and the Hankel function of the first kind#^ for polar co-ordinates. 
The even and odd angular Mathieu functions (Se m , So m ) are analogous to the cosine and sine functions, 


respectively, in polar co-ordinates. The normalization factors, an ^ ^ m are defined by Stratton (1941, 
equations (18) and (19) pg. 376). 

If the incident angle 4>o is set to 0, equation (5) simplifies to 
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This exact solution (eq. (5) or (6)) is used to validate the modal element solution. In general, the 
scattered field is represented by 



(7) 
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m = 0 


where the unknown coefficients A + m and are to be determined by solving the global matrix equation. If the 
incident angle <j>o = 0, then the odd terms vanish (2?+ m = 0). If the incident angle is nonzero, or if the scatterer 
is not symmetric about the x axis, then the odd terms must be included. 

Formulas to estimate the number of modes M coe( needed for convergence for scattering from circles 
can be found in Baumeister and Kreider (1993a). As the eccentricity of the ellipse increases, this number 
increases to a limit of roughly twice the value for circles. 

The numerical calculation of Mathieu functions is challenging. There are several sets of accepted 
notation and several types of normalization of the functions. The functions are represented as a series, each 
term of which depends on an eigenvalue that must be calculated by solving an equation involving a continued 
fraction. Details of the theory of Mathieu functions may be found in Magnus, et al., (1966). Details of the 
numerical calculation of Mathieu functions may be found in Blanch (1964), Canosa (1971), Clemm (1969), 
Hirsch (1972) and Toyama and Shogen (1984). Clemm’s algorithm, as improved by Hirsch, was adapted for 
use in this study. 

INTERFACE CONDITIONS 

At the interface S between the finite element domain and the homogeneous domain, both pressure and 
velocity are continuous (Temkin, 1981, pg. 80). The continuity of pressure 

“ p\s~ ~ ^ ^ 

is expressed numerically by an integral weighting procedure, so the continuity of pressure at the interface is 
represented by the set of equations 


0=2 n 

fs W m .[Pa-p] ds = 0 m* = 0,1.2,. ..,M coef -\ (9) 

e=o 

There must be one weight equation for each unknown in eqn (7); if odd terms are included then 
additional weights So m (c,r |) are introduced. 

The choice of weight functions is crucial, since, as discussed under Results, the solution is quite 
sensitive to the weights used. The natural choice is to use W m = Se m , since the angular Mathieu functions are 
orthogonal over the surface of the ellipse (Abramowitz and Stegun (1972), section 20.5). However, this choice 
does not give perfect results, so that possibly a better choice of weights could be found, although after some 
experimentation, it seems unlikely. 

The integration element ds in equation (9) is given by 


ds 2 - r 2 d& 2 + dr 2 
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where the radius at any point on the ellipse is 
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The discretization of these equations, as well as the imposition of the continuity of the pressure 
gradient across S , is done as in Baumeister and Kreider (1993a). 


FINITE ELEMENT SOLUTION 

The finite element domain is divided into triangular elements with unknown acoustic pressure at the 
nodes. It is assumed that all material properties are constant in each element. The setup of the global matrix 
follows the derivation in Baumeister and Kreider (1993a). 

The weak finite element formulation of equation (3) contains a surface integral with the normal 
derivative of p, which is equated to the normal derivative of p a . The normal derivative of the scattered field 
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requires some explanation. Notice that the notation for the arguments of the Mathieu functions has been 
streamlined, and uses (u,v) co-ordinates instead of (£,/n), to be consistent with the computer code. It is 
convenient to write the normal derivative in terms of rectangular co-ordinates 


dp 5 dp s a dp s . o 
= -7— cos p + sin p 


dn dx 


dy 


(13) 


where p is the angle between the normal to the surface and the horizontal. The code can provide the values of 
dp s /du and dp s ldv , so by writing 


dp s _ 3 p s 3 u + dp s 3v 
3jc 3 u 3jc 3v dx 


(14) 
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it remains only to determine du/dx , duldy, dv/dx and dvldy. This is done by differentiating the transformation 
equations (1) each with respect to x and y to obtain 


(0 


(A 

0 

-B 

o 

0 


0 

A 

0 

-B 

0 


B 

0 

A 

0 

W 


lo 

B 

0 

A) 


f'd i£\ 
dx 
3 u 
dy 
3v 

dx 

dv 


(15) 


where A = d sinh u cos v and B = d cosh u sin v. This 4 by 4 block matrix is easily inverted to give the desired 
quantities. Substituting them into equation (14) and then equation (13) gives the final form of the normal 
derivative 
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where the primes signify derivatives with respect to the arguments. 

RESULTS 


( 16 ) 


Three examples were chosen to illustrate the behavior of the modal element method when applied to 
scattering from ellipses. In Example 1, the scattered pressure in the near field is shown for a range of 
eccentricities. Example 2 shows the near and far scattered fields for a higher frequency. In Example 3, the 
finite element solution is compared to the modal solution at the surface of the ellipse. 

In each of the examples, the ellipse is defined with semi-major axis a = 1 (oriented along the x axis), 
and the semi-minor axis b varies between 0 and 1 to change the eccentricity. Near and far field polar plots 
(Examples 1 and 2) are taken on circles of the indicated radius to conform to usual far field results, while the 
surface plots of Example 3 are taken on the surface of the ellipse. The incident plane wave travels in the +x 
direction with a fixed frequency. 

In the examples, the ellipse is taken to be hard, which is simulated numerically through an impedance 
mismatch induced by setting 8 = 1 - 10 19 i and \i = 1 for each internal finite element. This feature allows 
greater flexibility in the numerical implementation of the method, because penetrable or coated rigid bodies 
may be studied with only slight modifications to the computer code, mainly in grid generation. For penetrable 
bodies, the internal grid is more extensive, while for coated bodies, several rings in the coating region may be 
needed. 

For examples with wave number k = 1, there are 264 grid points and 11 modal terms used, while for 
those with wave number k = 5, there are 384 grid points and 16 modal terms used. For all examples, the 
incident angle is set to 0 and the thickness of the ring (fig. 2) is set to 0.001. 

Example 1 — Varying The Eccentricity Of The Ellipse. — The incident plane wave has wave number 
k = 1 . Figures 4(a) to (d) show the modal solution for the scattered pressure field in the near field for a set of 
ellipses of increasing eccentricity. The exact solutions are represented by solid lines. As seen in figure 4, the 
errors increase slightly as b gets smaller (thinner ellipse). However, even for the most elongated scatterer, the 
method successfully calculates the field. For this reason, it is believed that the calculation of the far field 
pressure field from a realistic airfoil shape would be trustworthy. 

Example 2 — Near And Far Fields At Higher Frequency . — The incident plane wave has wave 
number k - 5 with semi -major axis a = 1 and semi-minor axis b = 0.5. Figures 5(a) to (d)show polar plots of 
the modal solution around circles of increasing radius. The exact solutions are represented by solid lines. It is 
clear that the solution improves farther from the scatterer, which is a consequence of the fact that the 
convergence properties of equation (7) improve as ^ (or u) increases. If only the far field is desired, a 
relatively coarse grid and fewer modal terms can be used. In the near field, more terms are required to obtain 
the same level of accuracy (Baumeister and Kreider 1993a, fig 4). Of course, as the frequency increases, a 
finer finite element grid and more modal terms must be used, leading to a larger global matrix and longer 
computation times. 

Example 3 — Comparison Of Finite Element And Modal Solutions. — The incident plane wave has 
wave number k - 1 with semi-major axis a = 1 and semi-minor axis b = 0.8. Figures 6(a) and (b) show the 
scattered pressure field at the surface of the ellipse based on the finite element (fig. 6(a)) and the modal 
(fig. 6(b)), solutions. The exact solutions are represented by solid lines. Previously, as has been shown in 
figures 4 and 5, the modal element solution has been shown to be very accurate in the far field. Now, as seen 
in figure 6(b), the modal element solution is quite good in the near field on the surface of the ellipse. 
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However, the finite element solution is clearly inadequate. For more elongated ellipses, the finite element 
solution is much worse, while the modal solution shows only a slight degradation in quality. 

The reason for this disparity is that the finite element solution is sensitive to the choice of weights 
used in the interface condition (9). The natural choice is the set of angular Mathieu functions Se m > but when 
this set (or variations such as cosine functions or phase shifted Mathieu or cosine functions) is used, the 
interface condition is satisfied globally but not locally. This means that although 

X Pi - Pn (*, ) = 0 0 

boundary nodes 1 a 1 


it is not true that 


X I p t — p n (jc- )| = 0. 

boundary nodes ' 1 a 1 ' 


(18) 


The finite element solution oscillates around the exact solution in such a way that the global 
condition is satisfied. Since the modal coefficients depend on the global condition, the modal solution is good 
even when the finite element solution is not. 

Attempts were made to strengthen the interface condition. First, different weights, including cosine 
functions and various phase shifted cosine functions, were used without success. Next, the condition itself was 
modified, by adding a tangential derivative: 



(19) 


But, by integrating the new term by parts, it is easy to see that the net effect is slight — for cosine weights, for 
example, the new weight is simply a phase shifted cosine. In all the trials, there was no improvement gained 
by changing the weights. In addition, changing the number of boundary nodes used in the quadrature scheme 
(in general, it is not necessary to use them all) had virtually no effect on the results. 

For scattering from hard bodies, this error in the finite element solution is irrelevant, since the modal 
solution accurately portrays the scattered field. For scattering from soft bodies, however, the internal field 
may be important, and the algorithm must be improved. There are two possible approaches to resolve the 
difficulty. 

First, the following iterative scheme could be implemented. After solving for the modal solution on 
the boundary, the internal field could be solved by using that solution as a Dirichlet boundary condition, and 
solving only for the finite element unknowns inside the body. 

Second, the condition could be made more stringent by giving up linearity. Possible forms are 


' boundary nodes 


(Pi-Pa( x i>f m0 ’ X 


boundary nodes 


p,-p a ( x i) l =0 - 


( 20 ) 


The resulting equations could be solved by a standard nonlinear system solver, but at greater computational 
expense. 

The numerical results can be summarized as follows: 

1. As the eccentricity of the ellipse and/or the incidient frequency increases, the problem requires a 
finer finite element grid and more modal terms. For the problems considered, the total number of unknowns 
ranged from 275 to 400, and the results shown are quite good. 

2. Far field results are better than near field results in all cases, because the modal series (eq. (7)) 
requires fewer terms to converge in the far field than in the near field. 

3. The finite element solution on the scattering body is often much poorer than the modal solution. 
This is due to sensitivity in the form of the interface condition. For hard scatterers, this is not a problem, while 
for soft scatterers, the algorithm can be modified to obtain the internal field. 
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CONCLUSION 


A study of acoustic scattering by ellipses using the modal element method has been presented. The 
method yields acceptable numerical results for scattering from hard ellipses of any eccentricity for a range of 
moderate frequencies. For soft ellipses, the method must be modified to obtain the internal field. The 
numerical cases presented are used as validation for the method, and indicate that when applied to scattering 
from realistic airfoil shapes, the modal element method should yield accurate results. In addition, the airfoil 
can be embedded in an ellipse rather than a circle, decreasing the size of the finite element grid, and hence 
reducing the computing requirements for the problem. 
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Figure 1 . — Finite element grid system. 
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Figure 2. — Finite element ring grid system for rigid bodies. 
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Exact analysis, ABS P 
Numerical solution, ABS P 
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Figure 4. — Near field polar plot of total acoustic pressure around a rigid ellipse subjected to plane wave impingement recon- 
structed from calculated modal coefficients at specified b (k = 1 , a = 1 , r = 1 .0). (a) b = 1 .0. (b) b = 0.8. (c) b = 0.5. (d) b = 0.2. 
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Ilipse subjected to plane wave impingement reconstructed from 
1, b = 0.5). (a) r = 1. (b) r= 3. (c) r = 5. (d) r = 10. 
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